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ABSTRACT 

We use homogeneous samples of radio-quiet Seyfert 1 galaxies and QSOs selected 
from the Sloan Digital Sky Survey to investigate the connection between the velocity 
shift and the equivalent width (EW) of the [OlII] A5007 emission line, and their cor- 
relations with physical parameters of active galactic nuclei (AGNs). We find a sig- 
nificant and negative correlation between the EW of the core component, EW(core), 
and the blueshift of either the core (the peak), the wing, or the total profile of [O III] 
emission; it is fairly strong for the blueshift of the total profile particularly. How- 
ever, both quantities (EW and velocity shift) generally have only weak, if any, corre- 
lations with fundamental AGN parameters such as the nuclear continuum luminosity 
at 5100 A (L5100), black hole mass (Mbh), and the Eddington ratio (L/Lndd); these cor- 
relations include the classical Baldwin effect of EW(core), an inverse Baldwin effect 
of EW(wing), and the relationship between velocity shifts and L/L^dd- Our findings 
suggest that both the large object-to-object variation in the strength of [O III] emission 
and the blueshift-EW(core) connection are not governed primarily by fundamental 
AGN parameters such as L5100, Mbh and L/L^dd- We propose that the ISM conditions 
of the host galaxies play a major role instead in the diversity of the [O III] properties 
in active galaxies. This suggests that the use of [OlII] A5007 luminosity as proxy of 
AGN luminosity does not depend strongly on the above-mentioned fundamental AGN 
parameters. 
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1. Introduction 

In terms of size, the narrow-line region (NLR) of active galactic nuclei (AGNs; including 
Seyfert galaxies and quasars) represents the link between the inner structure of the AGN — the 
dusty torus and the broad-line region (BLR) — and the interstellar medium in the host galaxy. The 
NLR bridges the region dominated by the gravitational field of the central massive black hole and 
the region dominated by the gravitational potential of the bulge stars of the host galaxy. To gain 
a complete understanding of the AGN phenomenon, the role of the NLR therefore needs to be 
understood. Unlike the BLR, the NLR, or at least the outer parts of the NLR, can be spatially 
resolved in nearby AGNs. Despite this, however, the NLR remains remarkably poorly understood 
and there are many unanswered questions about the structure and kinematics of the NLR and its 
relationship to the BLR/torus. 

[O III] A5007 is one of the most important NLR lines. Its equivalent width (EW) in the AGN 
ensemble varies dramatically by a factor of > 300 (from < 0.5 to 157 A; Baskin & Laor 2005b), 
one of the things showing the most object-to-object variation in AGNs. It is a dominant variable in 
the set of correlations making up the first principal component (commonly called "Eigenvector 1"; 
EVl) in the principal component analysis of QSO properties of Boroson & Green (1992); EVl is 
believed to be linked to certain fundamental parameters of the accretion process. [O III] luminosity 
is commonly used as a surrogate for the bolometric luminosity (e.g., Heckman et al. 2004) and its 
Doppler width as a surrogate for the stellar velocity dispersion (cr^) of the host galaxy (e.g.. Nelson 
& Whittle 1996, Wang & Lu 2001, Komossa & Xu 2007, Gaskell 2009a). 

Although the spatial distribution of [O III]-emitting gas in the NLR of AGNs is rather inho- 
mogeneous (see, e.g.. Das et al. 2006, Crenshaw et al. 2010), there are some global regularities in 
its strength and kinematics. Its EW was found to correlate negatively, albeit with a large scatter, 
with the continuum luminosity (Grindlay et al. 1980, Steiner 1981, Dietrich et al. 2002, Kovacevic 
et al. 2010; cf. Table 3 of Dong et al. 2011). This is similar to the famous "Baldwin effect" 
(Baldwin 1977, Baldwin et al. 1978) first discovered for the broad C IV A1549 line. The [O III] line 
profile can be regarded as consisting of two components: a narrow component (hereinafter, the 
'core') with a velocity close to the systemic redshift of the host galaxy and a low, broad base that is 
usually blueshifted (hereinafter, the 'wing'). There is mounting evidence that the core component 
is a good tracer of the gravitational potential of the host galaxy bulge (see, e.g., Gaskell 2009a and 
references therein) while the wing component arises from AGN outflows (see, e.g., Crenshaw et 
al. 2010, Komossa & Xu 2007 and references therein). Scattering could also be a contributor to 
broadening of the wings (Gaskell & Goosmann 2008). 

Besides these general global trends in the AGN ensemble, there are some interesting peculiar 
[OlII] phenomena that have been discovered in some particular AGN subclasses. For instance, in 
narrow-line Seyfert 1 s (NLS 1 s) that have broad H/3 line with full width at half maximum (FWHM) 
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< 2000 kms"^ (the "population A" of the Eigenvector 1 domain, Zamanov et al. [2002]; see Xu 
& Komossa [2010] for a recent review), [OlII] A5007 is generally very weak. Moveover, in a 
small fraction of NLSls, the whole [O III] A5007 profile has a bulk velocity shift of > 250 kms"^ 
blueward (the so-called "blue outliers", Zamanov et al. 2002, Boroson 2005, Zhou et al. 2006, 
Komossa et al. 2008) and there seems to be no non-blueshifted component from a canonical NLR 
as seen in low-ionization lines (Komossa et al. 2008). Komossa & Xu (2007) found that although 
the width of the core component of [O III] A5007 is generally a good surrogate for a^,, this however 
does not hold in "blue outliers". 

It is an interesting question whether or not these various regularities and phenomena can be 
unified and this has been looked at by previous authors. For instance, Zamanov et al. (2002) sug- 
gested that "blue outliers" with extreme blue wing components originate in strong AGN outflows, 
and moreover, are associated with the same outflows as the blue-shifted component of the broad 
Civ a 1549 emission line. The origin of the blueshifting of high-ionization broad emission lines 
such as Civ A 1549 remains a subject of debate. There is general agreement that there has to be a 
radial component of motion and some opacity source. Gaskell (1982) proposed that the blueshift- 
ing was due to an outflow of the high-ionization gas and obscuration by the accretion disk. While 
this has remained the most popular model, it has the major problem that outflows of the domi- 
nant line-emitting gas seem to be strongly ruled out by velocity-resolved reverberation mapping 
(GaskeU 1988, Koratkar & Gaskell 1989, see GaskeU 2009b for a review). GaskeU & Goosmann 
(2008) therefore argue that the blueshifting is due instead to scattering of photons by material with 
a net inward motion as suggested by Corbin (1990). 

Observationally, the C IV A 1549 Baldwin effect and blueshifting are related in the sense that 
higher-luminosity AGNs have both lower CiV EWs and larger blueshifts (Corbin 1991, 1992). 
This is further supported by the study of Richards et al. (2002). If this is also true for [O III] A5007, 
we will gain important insights into the origin and kinematics of the NLR. It might point to similar 
underlying causes for the two effects in both the BLR and NLR. Meanwhile, understanding these 
effects should enable us to calibrate better the Doppler widths and luminosities of [O III] A5007 to 
serve as proxies for the bulge gravitational potential and AGN bolometric luminosity. 

Motivated by these considerations, we present here the results of a systematic study of the 
velocity shift of [O III] A5007 emission line, its connection to the Baldwin effect, and connection 
to possible physical drivers. This study takes advantage of the unprecedented spectroscopic data 
from the Sloan Digital Sky Survey (SDSS; York et al. 2000). Throughout this paper, we use a 
cosmology with Hq = 70 km s"^ Mpc"\ = 0.3, and I^a = 0.7. 
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2. Sample and Measurements 
2.1. Samples 

Type 1 AGNs are chosen from the sample of Seyfert 1 galaxies and quasars of Dong et al. 
(2011) taken from the Fourth Data Release (Adelman-McCarthy et al. 2006) of the SDSS. We 
select type 1 AGNs having the [O III] A5007 emission line in the SDSS bandpass and with con- 
tinuum and emission lines suffering only minimally from contamination by host galaxy starlight. 
To ensure reliable analysis of the [O III] A5007 profile, we also require high spectral quality in the 
[O III] A5007 region. Our primary goal is to study the connection between velocity shifts and EWs. 
To ensure the reliable analysis of the velocity shift of [O III], we first construct a basic sample by 
taking the [SlI] AA6716,6731 doublet lines as the fiducial reference for the systematic redshift. 
[S II] is free from broad line contamination so using it is more accurate than using narrow H/3 (see 
Section 2.2). 

The final criteria used for the basic sample are therefore: (a) z < 0.3, to ensure that [S II] 
emission lines, as well as [O III] A5007 and broad H/3, are in the bandpass; (b) a median signal-to- 
noise ratio (S/N) > 10 per pixel for the optical spectra, and particularly > 15 in the [OlII] region 
(4995-5020 A); (c) the S/N of either [S II] A6716 or A6731 line > 3 to ensure [S II] to be a reliable 
reference for velocity shift; and (d) weak stellar absorption features, such that the rest-frame EWs 
of CaK (3934 A), CaH -i- He (3970 A), and U6 (4102 A) absorption features are undetected at 
< 2(7 significance. The last criterion ensures the measurement of the AGN luminosity and the 
emission-line EWs suffering little from the contamination of the host galaxy starlight (see the 
Appendix of Dong et al. 201 1). 

After removing duplications and sources with too many bad pixels in the H/3 -i- [O III] region, 
we obtain 565 type 1 AGNs. Radio jets may enhance [O III] emission (Labiano 2008) and thus 
possibly affect the results. We therefore use the method of Lu et al. (2007) to exclude sources 
detected by the FIRST (Faint Images of the Radio Sky at Twenty Centimeters) survey by matching 
with FIRST catalog (Becker et al. 1995). After this, our basic sample consist of 383 sources. 

Because of the redshift cutoff of z < 0.3 (to ensure the presence of [S II] as the redshift ref- 
erence), the dynamic ranges of AGN parameters is somewhat restricted (see Subsection 3.2). So 
in the investigation concerning [O III] EWs only, we further define an extended sample by relaxing 
both the redshift cutoff and the S/N requirement for [S II]. The extended sample comprises 1951 
radio-quiet type 1 AGNs at z < 0.8 culled from the Dong et al. (2011) sample, with a median S/N 
> 15 pixel"' in the [OlII] region. 
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2.2. Analysis of spectra 

Data analysis methods are as described in detail in Dong et al. (2011). We only provide a 
brief description here. Following Dong et al. (2008), we fit simultaneously the AGN featureless 
continuum, the Fell multiplets, and other emission lines using a code based on the MPFIT pack- 
age (Markwardt 2009), and if necessary we refine the emission-line fitting with the continuum 
subtracted spectra using a code described in detailed in Dong et al. (2005). The AGN continua are 
modeled by local power-laws for 4200-5600 A and the Ha regions. The Fe II emission is modeled 
with two separate sets of templates in analytical forms, one for the broad-line system and the other 
for the narrow-line system, using the identification and measurement of Fe II lines in I Zw 1 from 
Veron-Cetty et al. (2004), as listed in their Tables Al and A2. Within each system, we assume that 
the respective sets of Fe II lines have no relative velocity shifts and have the same relative strengths 
as those in I Zw 1 . Broad Fe II lines are assumed to have the same profile as broad H/?, while each 
narrow Fe II line is modeled with a Gaussian. During the fitting, the normalization and redshift of 
every system are set to be free parameters. Broad Hydrogen Balmer lines are fitted with as many 
Gaussians as is statistically justified. 

All narrow emission lines are fitted with a single Gaussian, except for the [O III] AA4959, 5007 
doublet lines. Each line of the doublet is modeled with two Gaussians, one accounting for the line 
core and the other for a possible blue wing, as seen in many objects. The Gaussian that has a smaller 
velocity shift relative to [S II], whether blueward or redward, is taken to be the core component. 
It also turns out to be narrower than the other Gaussian component (if the other one was present). 
The velocity shifts of the core component with respect to [S II] in the basic sample have a mean 
of -47 kms"^ (negative values denoting blueshifts) and a standard deviation of 72 kms"\ while 
the mean shift of the wing component is -225 kms"^ with a standard deviation of 240 kms"\ 
respectively. On average the core components comprise 54 per cent of the total emission, with 
a standard deviation of 0.16 dex. We list in Table 1 the fluxes, EWs and velocity shifts of the 
core and wing components and the whole profile of [O III] A5007 for all the objects in the basic 
sample. The velocity shifts are calculated by assuming that the [S II] doublet lines are located at 
the systematic redshift. In Table 2 we list the fluxes and EWs of the core and wing components of 
[O III] A5007 for all the objects at higher redshifts that are thus not in the basic sample but in the 
extended sample. The other spectral parameters of this sample are given in Dong et al. (2011)0 
An example is illustrated in Fig. 1 of fitting the [OlII] A5007 line profile of SDSS spectrum of 
J08013 1.58+354436.4 (i.e., the decomposition of the core and wing components). 

In the basic sample the velocity shifts of narrow H/3 with respect to [S II] have a mean of -21 
kms"' and a standard deviation of 82 kms"^ Such a scatter is a bit larger than that of the core 



'Also available online at |http://staff.ustc.edu.cn/~xbdong/Data_Release/ell_effect/ . 



-6- 



component of [OlII]. So if we use narrow H/3 as the fiducial reference for the systematic redshift 
to study the velocity shift of [O III], it would smear the correlations of interest in this study. This 
should be especially true for the sources at relatively higher redshifts in the extended sample. Those 
sources have lower spectral S/N and a weaker narrow H/3 component than the sources in the basic 
sample, thus the deblending of narrow H/3 from the broad component has a larger uncertainty. The 
use of narrow H/3 is checked in Section 3. 

We notice that there is usually a prominent Fe II emission feature immediately redward of 
[OIII]A5007, which is dominated by Fell 42 A5018 and [FeII]20F A5020. To ensure that the 
[O III] measurements suffer little from poor subtraction of this Fe II feature, we visually inspect the 
fits for the entire extended sample and fine-tune the fits for some objects by carefully matching 
the model to this local feature. In such objects, the relative strengths of Fe II lines are significantly 
different from those of I Zw 1 on which the Fe II template data are based. 



3. Results 

3.1. The connection between the blueshifting and Baldwin effects for [O III] A5007 

We explore the connection between the velocity shift and strength of [O III] A5007 in the 
basic sample. We calculate the velocity shifts of the core and wing components and the centroid 
(i.e., intensity-weighted average) of the whole profile ( Av(core), Av(wing), and Av(centroid), 
respectively) with respect to [SlI]. The line strength is characterized by the rest- frame EW, as 
measured by dividing the emission line flux by the underlying AGN continuum flux density at the 
central wavelength, for the core and wing components and the whole [O III] emission, respectively. 
We perform bivariate Spearman rank correlations to test the relationship between these kinds of 
velocity shifts and EWs. This method tests for not only a linear relation but a monotonic one. 
The results are summarized in Table 3, where we report the Spearman coefficient (r^) and the 
probability (Pnuii) that a correlation is not present. 

A striking regularity emerges from the correlation matrix that the velocity shifts of the core 
and wing components and the whole profile, correlate significantly with the EWofthe whole [O III] 
emission (Pnuw ^ 10"^ j and particularly with that of the core component, EW(core) (Pnuii < 10"^ j. 
The strongest correlation is between EW(core) and Av(centroid), the velocity shift of the whole 
profile (rs = 0.5 1 and Pnuii ^ 10"^^). For any a specific kind of velocity shift, the correlation is most 
significant with EW(core), yet not very significant (or not significant at all) with the EW(wing) (see 
Table 3). For the EWs of both the core component and the whole emission, the correlation is more 
significant with Av(centroid) rather than with Av(wing) or Av(core), yet their correlations with 
Av(core) are still very significant (Pnuw ^ 10"^). 
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In Fig. 2, we show the relationships between the EW of the core component and the velocity 
shifts of core component and the whole profile. Note that our Av(core) is effectively the same 
as the velocity shift measured from the peak of the whole [O III] profile ( Av(peak) ), because 
by fitting the core component one is always getting the peak of the whole [OllI] profile as the 
central wavelength of the 1 -Gaussian model. Av(peak) is commonly used in the definition of 
"blue outliers" in the literature (e.g., Zomanov et al. 2002, Boroson 2005, Komossa et al. 2008). 

To illustrate the [O III] velocity shift-EW relationship, we divide the entire sample into three 
subsamples according to Av(core) [namely, Av(peak)] with Av(core) < -150 kms"\ -100 ^ 
Av(core) < -50 kms"', and Av(core) ^ -50 kms"^ We then construct a composite spectrum 
for each subsample. We first subtract broad H/3 and Fell emission from the SDSS spectra, nor- 
malize every spectrum by the average flux density in the continuum window around 4200 A, and 
then construct the composite in the same manner as Richards et al. (2002) (see their Fig. 4). The 
composite spectra are demonstrated in Fig. 3 (panel a). 

The variations in EWs and velocity shifts among the composite spectra can readily be seen 
to be consistent with our correlation analysis in the sense that the more the peak of the line is 
blueshifted, the more the EWs of the total emission and, particularly, the core component, decrease 
dramatically, while the blue wing changes much less. This is similar to the case of C IV as Richards 
et al. (2002) found (see also Richards et al. 2010, Wang et al. 2011). In addition, there seems 
to be a trend of a positive connection between the blueshift of the peak and the strength of the 
wing component, which is not very significant according to the above correlation result between 
Av(core) and EW(wing) yet {P^„i\ = 4%). 



3.2. Correlations with physical AGN parameters 

It is interesting to explore whether the blueshifting/Baldwin effect of [O III] A5007 are driven 
by some fundamental physical parameter. We have therefore investigated the correlations of the 
velocity shifts and EWs of the core and wing components and of the total [O III] emission with 
broad H/? FWHM, continuum luminosity L5100 = ALa(5100 A), Mbh, and L/L^dd, for the sources 
in the basic sample. We calculate the black hole masses based on the FWHM of broad H/? using 
the formalism presented in Wang et al. (2009, their Eqn. 11). This formalism was calibrated with 
recently updated reverberation mapping-based masses and assuming BLR radius R oc L"^ (Bentz 
et al. 2009). The Eddington ratios are calculated assuming a bolometric luminosity correction 
^boi ~ 9 ALa(5100 A) (Elvis et al. 1994; Kaspi et al. 2000). The correlation results are summarized 
in Table 4. 

We find that there are some significant correlations (Pnuu ^ 10"^). The EW of the core com- 
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ponent has a negative correlation with the continuum luminosity (i.e., the classical Baldwin effect) 
and with the Eddington ratio (see also Dong et al. 2011). There is also a positive correlation of 
the blueshifting with the Eddington ratio, but none of these correlations are strong (rg < 0.4). In 
particular, no correlation of any EW or velocity shift with these physical AGN parameters is more 
significant (stronger) than the Av(centroid)-EW(core) relationship for the same sample. 

We also divide the 383 sources into three subsamples by L5100 and three subsamples by L/L^dd- 
We construct composite spectra for each subsample as described in Section 3.1. The three L5100 
subsamples have L5100 > lO"^"*-^, lO"*"^-^ > L5100 > 10"*^ and L5100 < 10"^^^ erg s"\ respectively, and 
the three L/Ledd subsamples have L/Lndd > 10-°\ lO'^-^ > L/Lgdd > 10"° ^ and L/Ledd < 10-° ^ 
The composite spectra are shown in panels b and c of Fig. 3. All of the already mentioned 
correlations can be seen. We also plot the relatively strong correlations of the velocity shifts of the 
peak (i.e., the core) and centroid of [O III] A5007 with Eddington ratio (Fig. 4). 

In summary, the results worthy of note in the correlation matrix are as follows: 

• The most significant correlation in the EW-related correlation matrix is the positive one 
between the intensity ratio of the wing to the core component and the continuum luminosity 
(Pnuii = 8 X 10"'° for the 383 sources), which hints a positive correlation — albeit a weak 
one to be sure — between the EW of the wing component and the continuum luminosity. 
This is confirmed by our analysis using the extended sample free from the redshift cutoff 
required by the presence of [S II] (see below), which gives a significant correlation between 
EW(wing) and L5100 with Pnuw = 6 x 10"'^. The significance is high because of the large 
number of sources (1951) but is still not strong (r^ ^ 0.2). The correlation of the EW of 
the core component with luminosity is slightly weaker than that with Eddington ratio in the 
basic sample (and much weaker in the extended sample — see below). Both are much more 
significant than those with broad-line FWHM or black hole mass. 

• The most significant correlations of all the velocity shifts of the core and wing components 
and of the total profile are with L/L^dd (Pnuii < 10~^), rather than with broad-line FWHM, 
L5100 or Mbh- These correlations are in the sense that the higher the L/L^dd the more the line 
is blueshifted. Moreover, all the three kinds of velocity shifts have no significant correlation 
with Mbh at all (Pnuii > 10%). Furthermore, all the objects with large blueshifts have a high 
Eddington ratio, but the reverse is not true (see Fig. 4; also Komossa et al. 2008). 

The dynamic ranges of L5100, Mbh and L/L^dd in the basic sample are 2, 1.5 and 1.5 dex 
respectively, which are not very large. It is thus possible that some weak but real correlations are 
obscured by the significant measurement errors (see Section 4.1 below; also Wang et al. 2009, 
Dong et al. 2011). To check the correlations concerning EWs that do not require [SlI] as the 
reference for systematic redshift, we use the extended sample as defined in Section 2.1. This 



-9- 



sample enlarges the dynamic ranges of L5100, Mbh and L/L^^^ by ~ 1, 1, and 0.5 dex respectively. 
The results are also summarized in Table 4 (bottom panel). The main conclusions remains almost 
unchanged with similar correlation strengths, except for the EW of the wing component. As above 
mentioned, EW(wing) now has a very significant (yet still not strong) correlation with L5100. It 
has a less significant correlation with Mbh, and almost no correlation with L/L^a still. Another 
notable change is that the significance of the correlation of EW(core) with L/L^a now becomes 
much greater than that with L5100, although the correlation is still not tight. 

To check for possible effects of the formalism we use to estimated the black hole masses, 
we re-examine the correlation analysis with Mbh calculated using several other commonly used 
formalisms (e.g., Greene & Ho 2005, Collin et al. 2006, Vestergaard & Peterson 2006, Shen et 
al. 2010). All the tests give similar results to those listed in Table 4 (see also Table 3 of Dong et 
al. 201 1). This is mainly because the dynamic range of Mbh covered by our samples is not very 
large (~ 2.5 dex for the extended sample, centered at ~ 10^ Mq), and in this range the various 
formalisms for calculating the black hole masses have only subtle differences from one another 
(Wang et al. 2009). 

Furthermore, we consider possible effect that our sample selection and spectra analysis may 
bring. In Section 2.1, we simply use the limiting radio flux of FIRST to minimize the effect of 
possible AGN jets. It is possible that many of the excluded radio sources are not radio loud but 
just nearby and bright. To check for this, we perform the correlation tests using the 565 sources 
before the radio exclusion. It turns out that, the inclusion of the radio sources only increases the 
correlation coefficients by a few per cent (rarely by 10 per cent). Particularly, all the correlations 
we concem in this paper, such as the blueshift-EW(core) relationship, stay nearly unchanged. This 
is mainly because the radio sources account for only 30 per cent of the 565 sources. 

We check the effect of using narrow H/3 as the fiducial reference for the systematic redshift. 
When it is applied to the basic sample, the strengths of the correlations concerning [O III] velocity 
shifts as listed in Tables 3 and 4 are all reduced significantly but the relative strengths among 
the correlations keep almost unchanged. When it is applied to the extended sample, the existing 
correlations get more statistically significant due to the larger sample size, yet the strengths (rg) 
get weaker; e.g., = 0.34 (Pnuii ^ 10"'°") for the correlation of Av(centroid) with EW(core), and 
Ts = 0.29 (Pnuii = lO"''') for its correlation with L/L^dd- Again, the relative strengths among the 
correlations keep almost unchanged. 

It is, obviously, rather artificial to decompose the [O III] profile into "core" and "wing" com- 
ponents by using simply two Gaussians. There is, for example, no theoretical justification that the 
wing component should have a Gaussian profile. And, decomposing the line profile into differ- 
ent components certainly depends on the line width at the SDSS spectral resolution. For objects 
with small line widths and small velocity shifts, the decomposition would be more difficult. We 



do several investigations to check possible decomposition effects. We find that the width of the 
decomposed core component, which is the most sensitive one to the spectral resolution, has only 
weak (rs < 0.1) correlations with EW(wing), EW(total), and the velocity shifts of the core and 
wing components and of the total emission, and no correlation (Pnun = 30%) with EW(core) at all. 
The width of the decomposed wing component has no strong correlations with the EWs or velocity 
shifts of the core component and of the total emission too, and almost no correlations with either 
the EW or velocity shift (Pnuii > 1%) of the wing component. These mean that the decomposition 
process would not bring a serious systematic effect into our results. 



4. Conclusion and Discussions 

4.1. [O III] properties are not governed preponderantly by AGN parameters? 

Despite the crude nature of the decomposition, the primary results about the velocity shifts 
and the EWs of [O III] A5007 (and its core and wing components) appear quite robust, according to 
both the correlation analysis and the composite spectrum analysis. Our main result is that, while 
there is a fairly strong correlation between the EW of the core component and the velocity shifts, 
there are no similarly strong correlations of the EWs or velocity shifts of either the cores, wings, or 
of the whole emission with black hole mass, Eddington ratio, or AGN luminosity. Thus a tentative 
conclusion is that neither the large object-to-object variation of the EW of [O III] A5007 nor its 
blueshift are readily attributable to some fundamental AGN parameters such as L5100, Mbh and 
L/LEdd (see also Baskin & Laor 2005b). Likewise, neither does the velocity shift-EW relationship 
seem to be primarily driven by some AGN parameter. The weakness of the correlations of [O III] 
properties with fundamental AGN parameters (including nuclear luminosity, Mbh and L/L^dd) 
have also been reported in the literature by using samples with much larger dynamic ranges (e.g., 
Dietrich et al. 2002, Netzer et al. 2004). 

This is in sharp contrast with the results for the EWs of broad Mg II A2800 doublet emission 
lines and of narrow optical Fe II emission lines, and the intensity ratios of ultraviolet and optical 
Fe II to Mg II, as recently discovered by Dong et al. (2009a, 20 11). They found that those EWs and 
intensity ratios correlate with L/L^dd strongly {r^ > 0.5 and Pnuii ^ 10"'^; see the above references), 
by using the same data set as the present study. Considering the effects of measurement errors are 
similar on the correlations of both those lines and [O III], we believe that the correlations of [O III] 
properties with L5100, Mbh, and L/L^dd are intrinsically weak. For those emission lines arising 
from the BLR and inner NLR (narrow Fe II lines, Dong et al. 2010), Dong et al. (2009a, 2009b, 
2011) proposed that the essential physical mechanism is that Eddington ratio regulates the global 
properties (particularly, the distribution of column density) of the clouds gravitationally bound in 
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the line-emitting regions. It appears now that the mechanism operating in the BLR and the inner 
NLR does not play a major role in the canonical NLR located far out in the host galaxies and free 
from the gravitational influence of the central supermassive black holes. A useful positive inference 
is that, at least in a statistical sense, the above-mentioned use of\0 III] A5007 as the proxy ofAGN 
luminosity does not depend seriously on L5100, Mbh and L/L^i^^. 

On the other hand, although the [OlII] properties do not depend as strongly on L5100, A/bh 
or L/Lgda as the above-mentioned emission lines arising from the inner NLR and the BLR, there 
are still some significant (albeit not strong) correlations between [O III] properties and those funda- 
mental AGN parameters. The EW of the core component has the most significant (anti-)correlation 
with L/Lndd, a less significant correlation with continuum luminosity, a much lower correlation 
with the broad-line FWHM, and almost no correlation with Mbh- For the wing component of 
[O III] the most significant correlation is with continuum luminosity. The behavior of the total 
emission results from the combination of the above two factors, and is mainly dominated by the 
core component. The magnitude of the blueshifting of the core, the wing, and the whole profile 
has the most significant correlation with L/Ledd, rather than with the other parameters. Because 
the SDSS spectroscopic survey is magnitude-limited and, moreover, because both Mbh and L/L^dd 
are constructed from broad-H/3 FWHM and L5100, there are apparent (likely not intrinsic) correla- 
tions among these four quantities. For instance, the Spearman correlation coefficients of L/L^dd 
with broad-H/3 FWHM, L5100, and Mbh for our basic sample of 383 sources are r^ = -0.82, 0.30, 
and -0.51, respectively. In light of the serious inter-dependence among these four quantities, the 
correlations of EW(wing) with broad-H/3 FWHM, Mbh, or L/L^dd are probably a secondary ef- 
fect of the stronger (thus presumably intrinsic) correlation with L5100, while the correlations of the 
velocity shifts and EW(core) with broad-H/3 FWHM, L5100, or Mbh are a secondary effect of that 
with L/Lfidd- This is probable given that there are additional systematic uncertainties plaguing the 
estimated values of L/L^dd- One effect is the large uncertainties in virial BH masses, which can 
be a factor of 4 (la) statistically, and perhaps as large as an order of magnitude for individual 
estimates (Vestergaard & Peterson 2006; Wang et al. 2009). Another uncertainty comes from the 
bolometric correction assumed for L5100, which is definitely an oversimplification given the diverse 
spectral energy distributions of AGNs (Ho 2008; Vasudevan & Fabian 2009; Grupe et al. 2010). 
Yet, since the above bivariate correlations of [O III] properties with L5100, Mbh and L/L^dd are not 
strong, being weaker than the Av(centroid)-EW(core) relationship and much weaker even than 
the inter-dependence among the AGN parameters, we refrain from an over-discussion on this is- 
sue at present. The weakness of those correlations either reflects true scatters that are caused by 
multiple (and somehow independent) processes and are thus irreducible, or suggests that there are 
other dominant variable(s) yet to be identified (see the discussions on various possible physical 
variables/processes in Baskin & Laor 2005b, Risaliti et al. 2011 and Caccianiga & Severgnini 
2011). 
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4.2. The velocity shift-EW(core) connection 

The fairly tight anticorrelation between the blueshifts and the EW of the core component ap- 
pears more notable than the weak correlations with fundamental AGN properties. Interestingly, as 
analyzed in Subsection 3.1, this pattern is just like the situation for the broad C IV A 1549 emission 
line (Richards et al. 2002, Bachev et al. 2004, Richards et al. 2010), which has already been 
noted by some researchers (e.g., Zamanov et al. 2002). This suggests that both the blueshifting 
phenomenon and the classical Baldwin effect are tightly linked (Corbin 1991, 1992). Note that 
the classical Baldwin effect only concerns the core component. The wing component can even 
show an inverse Baldwin effect — the higher the AGN luminosity, the larger the EW of the wing 
component. This is similar to the situation of C IV A1549 (Richards et al. 2010, Wang et al. 201 1). 

In light of the remarkable similarity in the velocity shift-EW connection between [O III] and 
C IV, Occam's razor would suggest that the mechanism producing the two correlations is the same. 
A wind has been proposed as the cause of the C IV blueshifting (Gaskell 1982) and the correlations 
discussed above seem to be consistent with a scenario where the core of [O III] comes from the 
canonical extended NLR where the gas motions are dominated by the gravity of the bulge stars 
while the wing component arises from an outflow. This has been discussed by many authors (see 
Zamanov et al. 2002, Baskin & Laor 2005a, Richards et al. 2010, and Wang et al. 201 1 for details). 
In this model the "blue outliers" could just be extremes; presumably a higher relative accretion rate 
drives a stronger outflow. 

As pointed out in Section 4.1, the large object-to-object variation in the strength of [OlII] 
emission and the blueshift-EW(core) connection are not governed primarily by any fundamental 
AGN parameters such as L5100, Mbh and L/L^dd- It is likely that it is the ISM conditions of the host 
galaxies that mainly determine the diversity of the [O III] properties in active galaxies (particularly 
its strength; see Netzer et al. 2004). In the outflow scenario, it is natural to assume that the ISM 
gives rise to [O III] emission and decelerates the outflow. Although AGN activity determines the 
launching speed of outflow, the final speed of the outflow depends on the density and density 
profile of the ISM. Meanwhile, the density and the density profile also constrain the amount of 
[O III] -emitting gas and its emissivity. At the same AGN activity level, denser ISM will produce 
more [O III] emission while decelerate the outflow efficiently, leading to higher EW and lower 
velocity. On the other hand, the similar blueshift-EW relationship in C IV may be linked to the 
shape of the ionizing continuum (Richards et al. 2010). So it is possible that the above two factors 
work together to produce the blueshift-EW(core) connection. 

Yet, we should also mention that, while the outflow scenario is attractive, there seem to be 
problems. For instance, the outflow explanation for the blueshifting of C IV (Gaskell 1982) clashes 
head on with the BLR kinematics implied by velocity-resolved reverberation mapping. Velocity- 
resolved kinematics have long favored gravitationally-dominated motions (Gaskell 1988, Koratkar 
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& Gaskell 1989, 1991; Gaskell 2009b, 2010b) and apparent cases of outflow are unlikely to be real 
but are instead probably the result of off-axis variability (Gaskell 2010a). Gaskell & Goosmann 
(2008) have proposed instead that the blueshifting is the result of scattering off material with a 
slight net inflow. This naturally explains the correlation of blueshifting with accretion rate. If the 
NLR and BLR blueshiftings have a similar cause, the blueshifting of [O III] too might be caused 
by scattering, as has been proposed by Gaskell & Goosmann (2008). 

Finally, we note that while all the AGNs with large blueshifts have a high Eddington ratio, the 
reverse is not true (see also Komossa et al. 2008). Again, this is very similar to the situation for 
broad CiV A1549 (Baskin & Laor 2005a). Here we must stress once again (see Section 4.1) that 
the correlations between velocity shift and any fundamental AGN parameter such as nuclear lumi- 
nosity, black hole mass, and the Eddington ratio, are not strong (see also Boroson 2005, Komossa 
et al. 2008). Thus it is not clear if the velocity shift-EW(core) connection for [O III] A5007 can 
be reduced to the effect of a specific underlying physical process of the AGN activity. Anyway, 
whatever the cause(s) of the blueshifting of high-ionization NLR and BLR lines, the correlations 
discussed here provide further support for an intimate connection between the NLR and the BLR. 
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Fig. 1. — Results of the emission-line profile fitting applied to the H/3-[0lII] region of SDSS 
J08013 1.58+354436.4. We plot the continuum and Fell subtracted spectrum (black), the sum of 
all the best-fitting components (red), the fitted narrow (green) and broad H/3 (cyan), and particularly 
the decomposed core (pink) and wing (blue) components of [O III] A5007 line. The green dashed 
vertical line denotes the wavelength for [O III] A5007 as inferred from the [S II] line. 
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Fig. 2. — Distributions of velocity shifts of the core component (top) and the centroid of the whole 
profile (bottom) of [O III] A5007 versus the equivalent width of the core component, for the basic 
sample of 383 objects. 
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Fig. 3. — Panel (a): The composite spectrum of different Av(core) bins. The blue, green and red 
lines are spectra with decreasing Av(core). Panel (b): The composite spectrum of different L5100 
bins. The blue, green and red lines are spectra with decreasing Ljioo- Panel (c): The composite 
spectrum of different L/LEdd bins. The blue, green and red lines are spectra with decreasing 

J^/J^Edd- 
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Fig. 4. — Distributions of velocity shifts of the core component (top) and the centroid of the whole 
profile (bottom) of [O III] A5007 versus the Eddington ratio, for the basic sample of 383 objects. 
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Table 1 . Parameters of the Core and Wing Components of [O III] A5007 Emission Line: Basic 

Sample 



SDSS Name 


logF(core) 


EW(core) 


logF(wmg) 


EW(wing) 


Av(core) 


Av(wing) 


Ai'(centroid) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J0()0410.80-104527.2 


-15.02 


4.39 


-15.00 


4.64 


-77.52 


-608.9 


-351.3 


J0()0904.54-103428.6 


-14.24 


33.43 


-14.59 


14.93 


-1.431 


-88.4 


-28.3 


J001247.93-084700.4 


-14.75 


6.55 


-14.46 


12.91 


-25.71 


-186.9 


-132.7 


J001340.73-1 11 100.6 


-14.56 


15.29 


-15.13 


4.10 


-319.6 


-1134.4 


-492.9 


J001416.92+145038.4 


-14.94 


9.24 


-14.85 


11.26 


-130.5 


-483.8 


-324.9 


J003723.49+000812.5 


-14.82 


13.12 


-14.97 


9.35 


-37.27 


-422.4 


-197.9 


J004930.90+153216.3 


-14.08 


32.48 


-14.62 


9.24 


16.22 


-480.5 


-94.2 


J005118.27+135448.0 


-14.07 


26.38 


-15.26 


1.72 


2.772 


265.9 


18.8 


J005328.80-085755.0 


-14.69 


10.91 


-14.98 


5.51 


-32.51 


-250.1 


-105.7 


J005709.93+144610.2 


-14.02 


6.72 


-14.19 


4.57 


106.7 


-487.3 


-134.7 



Note. — Parameters for the 383 objects in the basic sample. Column (1) official SDSS name; 
Column (2) Flux of the core component; Column (3) Equivalent Width of the core component; 
Column (4) Flux of the wing component; Column (5) Equivalent width of the wing component; 
Column (6) Velocity shift of the core component, with respect to [S II] AA67 16, 673 1 doublet lines; 
Column (7) Velocity shift of the wing component, with respect to [S II]; Column (8) Velocity shift 
of the centroid, with respect to [S II]. (This table is available in its entirety in a machine-readable 
form in the online journal. A portion is shown here for guidance regarding its form and content.) 
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Table 2. EW and Flux of the Core and Wing Components of [O III] A5007 Emission Line: 

Extended Sample 



SDSS Name 


logF(core) 


EW(core) 


logF(wmg) 


EW(wing) 


(1) 


(2) 


(3) 


(4) 


(5) 


J000043.95-09 1134.9 


-15.27 


4.38 


-15.42 


3.11 


J000102.19-102326.8 


-14.39 


35.43 


-14.66 


18.9 


JOOOl 10.96-105247.4 


-14.99 


5.52 


-14.70 


10.8 


J000545.61+153833.8 


-14.88 


7.73 


-14.59 


15.0 


J000945.46+001337.1 


-15.73 


2.77 


-15.03 


13.7 


J001030.58+010006.0 


-13.84 


82.26 


-14.43 


21.4 


J001327.31+OO5232.0 


-14.33 


24.42 


-14.34 


23.6 


J001630.43-093853.4 


-14.98 


7.48 


-15.00 


7.18 


J001725.35+141 132.5 


-14.74 


5.28 


-14.52 


8.82 


JOO 1855. 22-09 135 1.1 


-15.08 


3.85 


-14.70 


9.39 



Note. Parameters for the 1568 objects that are not in the basic sample but in the extended sample. 
Column (1) official SDSS name; Column (2) Flux of the core component; Column (3) Equivalent 
Width of the core component; Column (4) Flux of the wing component; Column (5) Equivalent 
width of the wing component. (This table is available in its entirety in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding its form and content.) 



Table 3: Correlation Results between Velocity Shift and Equivalent Width ^ 





EW(core) 


EW(wing) 


EW(total) 


Av(core) 

Av(wing) 

Av(centroid) 


0.23 (4e-06) 
0.41 (6e-17) 
0.51 (2e-26) 


0.10 (4e-02) 
0.19 (2e-04) 
0.11 (3e-02) 


0.20 (7e-05) 
0.38 (le-14) 
0.41 (8e-17) 



" For each entry, we list the Spearman rank correlation coefficient (rs) and the probability of the 
null hypothesis that the correlation is not present (Pnuu) in parenthesis, for the 383 objects in the 
basic sample. 
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Table 4: Correlation Results between [O III] Quantities and Physical Parameters " 







FWHM(H/3'^) 


-^-5100 ^ 






Av(core) 


383 


0.28 (3e-08) 


-0.12 (le-02) 


0.17 (7e-04) 


-0.35 (3e-12) 


Av(wing) 


383 


0.17 (7e-04) 


-0.25 (le-06) 


0.03 (5e-01) 


-0.31 (le-09) 


Av(centroid) 


383 


0.22 (2e-05) 


-0.27 (6e-08) 


0.06 (2e-01) 


-0.36 (2e-13) 


EW(core) 


383 


0.14(5e-03) 


-0.26 (2e-07) 


-0.01 (9e-01) 


-0.28 (4e-08) 


EW(wing) 


383 


0.10(6e-02) 


0.02 (7e-01) 


0.08 (le-01) 


-0.07 (le-01) 


EW(total) 


383 


0.14(5e-03) 


-0.18 (3e-04) 


0.03 (6e-01) 


-0.23 (5e-06) 


EW(wing) 
EW(core) 


383 


-0.06 (3e-01) 


0.31 (8e-10) 


0.09 (8e-02) 


0.22(le-05) 


EW(core) 


1951 


0.11 (le-06) 


-0.16(3e-12) 


-0.02 (4e-01) 


-0.22 (6e-23) 


EW(wing) 


1951 


0.08 (5e-04) 


0.19 (6e-18) 


0.16 (9e-13) 


0.07 (le-03) 


EW(total) 


1951 


0.12 (7e-08) 


0.01 (6e-01) 


0.08 (3e-04) 


-0.11 (3e-06) 


EW(wing) 
EW(core) 


1951 


-0.05 (2e-02) 


0.31 (6e-44) 


0.14 (5e-10) 


0.28 (3e-36) 



" For each entry, we list the Spearman rank correlation coefficient (rg) and the probability of the 
null hypothesis that the correlation is not present (Pnuii) in parenthesis. 
^ The number of the sources used in every bivariate correlation tests. 

L5100 ^ ALa(5100A); the BH masses are calculated using the formalism presented in Wang et 
al. (2009, their Eqn. 11); Eddington ratios (L/Lndd) are calculated assuming that the bolometric 
luminosity Lboi ^ 9L5100. 



